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Abstract. We revisit integrable discretizations for the nonlinear Schrodinger 
equation due to Ablowitz and Ladik. We demonstrate how their main drawback, 
the non-locality, can be overcome. Namely, we factorize the non-local difference 
scheme into the product of local ones. This must improve the performance of 
the scheme in the numerical computations dramatically. Using the equivalence of 
the Ablowitz-Ladik and the relativistic Toda hierarchies, we find the interpolating 
Hamiltonians for the local schemes and show how to solve them in terms of matrix 
factorizations. 



1 Introduction 



In [jl j -0] the author pushed forward a new method of finding integrable 
discretizations for integrable differential equations, based on the notion of r— 
matrix hierarchies and the related mathematical apparatus. The main idea of 
this approach is to seak for integrable discretizations in the same hierarchies 
where their continuous counterparts live. 

In fact, this method is not quite new. The first integrable discretizations 
which can be treated as an application of this method go back as far as to 
the year 1976, to the work of Ablowitz and Ladik ||. In the present note 
we revisit the Ablowitz-Ladik discretizations, improving them both from the 
aesthetical (theoretical) and the practical (computational) point of view. 

In H Ablowitz and Ladik proposed the following remarkable system of 
ordinary differential equations: 

qk = qk+i - 2qk + qk-i - qkn(qk+i + qk-i) 

(1.1) 

h = -r k+1 + 2r k -r k - 1 + q k r k (r k+1 + r k - 1 ). 

It may be considered either on an infinite lattice (k G z) under the boundary 
conditions of a rapid decay {\q k \, ~~ ► as k — > ±oo), or on a finite lattice 
(1 < k < N) under the periodic boundary conditions (q = q^, r = r N , 
qN+i = qi, Tiv+i = Ti). In any case we shall denote by q (r) the (infinite- or 
finite-dimensional) vector with the components q k (resp. r k ). 

In |5|] the system ( |1.1| ) appeared as a space discretization of the following 
system of partial differential equations: 

qt = q xx ~ 2q 2 r, r t = -r xx + 2qr 2 (1.2) 

(to perform the corresponding continuous limit, one has first to rescale in 
(|1.1| ) t i — > e~H, q k t— > eq k , r k \— > er k , and then to send e — > 0). 

It is important to notice that upon the change of the independent variable 
1 1— > it, i — V— 1, the system (|1.2|) allows a reduction 

r = ±q*, (1.3) 

leading to the nonlinear Schrodinger equation 

- iqt = q xx =f 2\q\ 2 q- (1-4) 
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(In ( |1.3| ) and below we use the asterisque * to denote the complex conjuga- 
tion). The same reduction is admissible also by the Ablowitz-Ladik system 
fllTTD, leading to 

- iqk = qk+i - 2q k + q k -\ =F |<?fc| 2 (<?fc+i + qk-i)- (1-5) 

Ablowitz and Ladik found also a commutation representation for the sys- 
tem (|1 . 1| ) - a semi-discrete version of the zero-curvature representation: 

L k = M k+1 L k - L k M k (1.6) 

with 2x2 matrices L k , M k depending on the variables q, r and on the 
additional (spectral) parameter A: 

A q k 

L fc = L fc (q,r)=| |, (1.7) 

r k A- 1 



M k = M k {q, r) 



/ A 2 - 1 - q k r k -i Mk ~ A 1 q k -i 



;i.s) 



\ Ar fc _i - A V fc -A 2 + 1 + q k -ir k 
Note that the linear problem associated with the matrix Lk, 

^ fc+1 = L k V k , (1.9) 

is a discretization of the linear Zakharov-Shabat problem, associated with 
the system ( |1.2|) , 

(K q \ 

® x =\ (1-10) 
\ r -i( J 

In H Ablowitz and Ladik made also the next step in discretizing the 
system (|1.2j ): they constructed a family of time discretizations of the system 
( |1.1| ). Although it was not stressed very explicitly, this time their approach 
to discretization was fundamentally different: they did not modify the linear 
problem fll.Sf) any more, restricting themselves with a choice of a suitable 
(discrete) -time evolution of the wave function ty k . Hence the basic feature 
of the time-discretizations in is following: they admit a discrete analog 
of the zero-curvature representation, 

L k V k = V k+1 L k (1.11) 
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with the same matrix L k as the underlying continuous time system. (In 



1.11|) and below we use the tilde to denote the h— shift in the discrete time 
hz). In a more modern language, the maps generated by the discretizations 
in H belong to the same integrable hierarchy as the continuous time system 



1.1 




;i.l2) 



The results of may be formulated as follows. 

Proposition 0. Let the matrix L k be given by ( |1.7|) , and let the entries 
of the matrix 

I A k 

V k = 

V c k 

have the following X- dependence: 

A k = 1 + hX 2 A k 2) + hA k 0) + h\- 2 A { k 2 \ 

v k = i + hx 2 v {2) + hvf ] + h\- 2 v[- 2 \ 
B k = hXB { k 1} + h\- 1 B k r 1 \ 
C k = hXC { k ] + h\- x C { k ~ l) . 

Then the discrete zero-curvature equation ( |1.11| ) implies the following expres- 
sions: 

A k = l-ha + ha + (\ 2 -A k ) + h(a^\- 2 -5_q k r k _ 1 )A k , (1.13) 
V k = l + h6 -hS+{\- 2 -D k )-h{5-\ 2 - a-q k - x r k )k k , (1.14) 
B k = h(a + Xq k - 5 + X' 1 q k ^i) + h(5^Xq k - a_A~V-i)Afc, (1.15) 
C k = h(a + Xr k -i - 5+A~V) + h(8-\r k -i - «_A _1 f fc )A fc , (1.16) 

together with the equations of motion: 

(qk-Qk)/h = a + q k+ i - a q k - 5 q k + <5 + g fc -i - (a + q k A k+1 + S + q k D k ) 
+ (5-q k+1 + <x_g fc _i)(l - q k r k )X k 

(1.17) 

{r k -r k )/h = -S + r k+ i + 5 r k + a r k - a + f fc _i + (5+r k D k+1 + a+r k A k ) 
-(a-Ffe+i + 5_r fe _i)(l - q k r k )K k 
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Here ocq, a+, a~, So, d~+, 5- are constants, and the functions A k , D k , A k 
satisfy the following difference relations 

A k+1 - A k = q k+ ir k - q k r k -i, (1-18) 

Dk+i ~ D k = q k r k+1 - q k - X r k , (1.19) 

Ajfe +1 (l - q k r k ) = A fc (l - q k r k ). (1.20) 

Remark. In the case of the rapidly decreasing boundary conditions 
there exists a canonical way to single out certain solutions of the difference 
equations ( |1.18| )-( |I.20| ) above, namely, by the conditions 



At, D k — ► 0, As. — > 1 as k — > ±oo, 



which results in 



fc-i 

A k = q k r k _ x + (Wj-i ~ (!-21) 

j=-oo 
k-l 

D k = q k „ir k + (?i-i r i _ I 1 - 22 ) 

j=-oo 

A, = ff (!- 23 ) 

In the case of the periodic boundary conditions to choose a particular solution 
one can use the same formulas with the sums and product starting from j = 
instead of j = — oo. 

The numbers a , 8 , a_, <5_, playing the role of constants of summation, 
are defined as soon as certain solutions A k , D k , A k have been fixed. 

The formulas ( |1 . 1 7| ) , ( |1.21| )-( |1.23| ) define the map which we shall denote 

T Ah (h; a , a+, a_; 8 , 6+, <L) : (q,r) h-> (q,r) 

It is important to notice that for the pure imaginary values of h this map 
obviously allows the reduction 

r = ±q*, 
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provided A k , D k , A k are chosen as in ( |1.21| )-( |I.23| ), and 



So = q*q, 8 



a* + , 5- 



;i.24) 



so that that there remains a three-parameter family of difference schemes 
satisfying this condition. 

The expressions ( 1.21 )— ( 1.23 ) serve as a source of a non-locality of the 
difference scheme, which is its major drawback. This feature makes any 
numerical realization of the numerical scheme extremely time-consuming. 

The numerical experiments reported in showed that even despite this 
drawback the Ablowitz-Ladik difference schemes are the best among the 
class of finite difference methods, being surpassed only by certain spectral 
numerical methods. 

In the present note we shall demonstrate how to factorize the non-local 
scheme ( |1.17|) into the product of very simple (in particular, local) schemes, 
which surely can speed up the performance of this scheme considerably. 



2 Ablowitz— Ladik hierarchy 
and its simplest flows 

From the modern point of view, the Ablowitz-Ladik system ( |1.1[ ) is a repre- 
sentative of a whole hierarchy of commuting Hamiltonian flows. Considering, 
for notational simplicity, the finite dimensional case, we define the Poisson 
bracket on the space R 2Ar (q, r) by the formula 

{%>rj} = (1 - q k r k )5 jk , {q k ,Qj} = {r k ,rj} = 0. (2.1) 

The Hamiltonians of the commuting flows are the coefficients in the Laurent 
expansion of the trace tr T/v(q, r, A) where T/v is the monodromy matrix 

Tjv = Ln ■ Ln-i ■ . . . ■ L 2 ■ Li, (2.2) 

supplied by the function 

N 

#o(q,r) = logdetTjv = £ log(l - q k r k ). (2.3) 

k =i 
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The involutivity of all integrals of motion follows from the fundamental r- 
matrix relation: 



{L(A)®L( / i)}=[L(A)®L( / i),p(A,/i)] ; 



(2.4) 



where 



/ 1 A 2 + M - 
' 2 A 2 — /x ; 




v o 



1 

2 

A 2 -/t 2 







_ 1 

2 









1 A 2 +^ 2 

2 A 2 - M 2 



(2.5) 



It is easy to see that the following two functions belong to the involutive 
family generated by tr T^: 



N N 

H+((l, r) = Qk+irk, H^(q, r) = ^ qkn+i- 

k=l k=l 



(2.6) 



The corresponding Hamiltonian flows are described by the differential 
equations 



^4 



Qk 



Qk 



q k+1 (l - q k r k ), f k = -r k -i(l-q k r k ), 
q k -i(l - q k r k ), r k = -r k +i(l - q k r k ), 



(2.7) 
(2.8) 



The flow generated by the Hamiltonian function ( |2.3| ) is described, up to the 
factor 2, by the differential equations 



Jo : <?fc = -2q k , r k = 2r k 



(2.9) 



The Ablowitz-Ladik flow proper is an obvious linear combination of these 
more fundamental and simple flows. According to the general theory M, each 
of the flows F±, Tq is described by the zero-curvature representation ( |1.11| ) 
with the same matrix L k , but with different matrices M k . The corresponding 
matrices M k are given by: 



^4 



M; 



(+) 




(2.10) 
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(-) 



M, 



( -A-V-i 

\ -A _1 r fc -A -2 + q k -\r k 
( -1 \ 



(0) 



(2-11) 



(2.12) 



3 Local discretizations for T- 



± 



We demonstrate now an unexpected fact. Namely, let only one of the four 



parameters a+, <5+, a_, <5_ of the Ablowitz-Ladik scheme not vanish (so 
that the resulting scheme approximates one of the flows T± rather than the 
original system (|1.1|) ). Then it is possible to render the scheme local. This 
results in four integrable maps (q, r) i— ► (q, r) which are described by local 
equations of motion and belong to the Ablowitz-Ladik hierachy, i.e. admit 



commutation representations with the matrix L k from ( |1.7|) . The commu- 
tation representations for the maps given in the following four Propositions 
could be proved by an easy and direct check, but we prefere to trace back 
the relations between our maps and the original formulation by Ablowitz and 
Ladik. 

Proposition 1. The Ablowitz-Ladik scheme ( |1.17 ) with the parameters 



a = 5q 



OC- 



5_ = 5 A 



0. 



OC- 



1 



is equivalent to the following map: 

- <lk)/h = gfc+i(l - qkh), 

(ffc - r k )/h = -f fc _i(l - q k r k ) 
approximating the flow T + . This map has the commutation representation 



T + {h) : 



(3.1) 



T + (h) : L k vt +) = V£{L 



with the matrix 



v k i+) = vr(^h) 



c+)< 



l + h\ 2 - hqkh-i h\q k 



h\r k . 



(3.2) 
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Proof. According to the Proposition 0, the matrix V k + ^ for the scheme 
T Ah (h; 0,1,0; 0,0,0) has the form 

1 + h\ 2 — hA k h\q k 

h\r k -t 1 

while the equations of motion read: 

(q~k - <lk)/h = q k +i - qkA k+1 , (r k - r k )/h = -r k -i + r k A k . (3.3) 

Here A k is the solution of the difference relation 

A k+ i - A k = q k+ ir k - q k r k -i (3.4) 

tending to as k — > ±oo in the case of the rapidly decaying boundary 
conditions. The Proposition will be demonstrated if we prove the following 
formula for A k . 

A k = q k r k -i (3.5) 

which comes on the place of the non-local expression ( |1.21[ ). To this end 
multiply the first equation in Q3.3|) by r k ^i, the second one by q k+ i, and add 
the two resulting equations: 

q k r k -i - qk+\r k + q k+ ih - q k r k -i = hq k+1 r k A k - hq k r k -iA k +i. 

Using fl3~4f ) we obtain: 

-A k+1 + A k + q k+1 r k - q k r k -i = hq k+1 r k A k - hq k r k -iA k+ i, 

which is equivalent to 

1 - hA k 1 - hA k+ i 



1 - hq k r k -i 1 - qk+ir k 



const. 



As the both quantities A k , q k r k _i tend to by k — > ±oo, this constant has 
to be equal to 1, which ends the proof. 

Proposition 2. The Ablowitz-Ladik scheme ( |1.17| ) with the parameters 

«o — <5o — a - — — a + — 0, 8 + = 1 
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is equivalent to the following map: 

!(<}k -Qk)/h = g r fc _i(l - q k r k ), 
(3.6) 
(h - r k )/h = -r fc+ i(l - q k r k ) 

approximating the flow This map has the commutation representation 



T-(h): L k vt ] = V^>L k 
( 



u>/t/i i/ie matrix 

<-) - T/H/ 



1 -/iA 1 g fc _ 1 



\ -/iA Vfe 1 — /iA 2 + hqk-iTk 



(3.7) 



Proof of this Proposition is completely parallel to that of the previous 
one and is therefore omitted. 

Proposition 3. The Ablowitz-Ladik scheme ( 1.17|) with the parameters 

ctQ = So = a + = 5 + = 5- = 0, a- — 1 

is equivalent to the following map: 

(Qk - Qk)/h = q k -i(l - q k r k ), 



(3.8) 



(r k - r k )/h = -r fe+ i(l - q k r k ) 
approximating the flow T-. This map has the commutation representation 



T~\-h) : L k wt ] = Wt+[L k 



with the matrix 



1 



/ 1 + hX 2 - hq k _ x r k -hX l q k -i 



1 - hq k -ir k 



\ 



-hX 1 r k 



(3.9) 

Proof. According to the Proposition 0, the matrix wj: ^ for the difference 
scheme 7al(^; 0; 0; 1> 0; 0; 0) h as the form 



W, 



(-) 



1 + hX~ 2 A k -hX- l q k ^A k 
-/iA _1 r fc A fc 1 + /ig fc _if fc A fc 
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The equations of motion for this scheme may be presented as 

(Qk~Qk)/h = g*-i(l-gfcffc)A fc , (f h -r k )/h = -f fc+1 (l-g fc r fc )A fe+1 . (3.10) 

Here is the solution of the difference equation 

A fc+ i(l - q k r k ) = A fc (l - q k r k ) (3.11) 

tending to 1 by k — > ±oo in the case of the rapidly decaying boundary 
conditions. It is easy to check that the proposition will be proved if we 
demonstrate the following formula for A k : 

A k = - * = (3.12) 
1 - hq k _ x r k 

which replaces in this case the general non-local expression (|1.23| ). To do 
this, we first re-write ( |3 . 1 U| ) as 

q k (l + hq k - X f k A k ) = q k + hq k - X A k , r k (l + hq k r k+1 A k+ i) = r k + hr k+1 A k+1 . 

Now multiply the first of these equations by r k , the second one by q k and 
subtract the two resulting equations: 

(1 - q k r k ){\ + /ig fc _if fc A fc ) = (1 - q k r k )(l + hq k r k+1 A k+1 ). 

According to ( p.ll| ), this is equivalent to 

Afc+i _ 1 + hq k r k+1 A k+1 
A k 1 + /ig fc _if fc A fc 

° r I I 

— + hq k -ir k = h hq k r k+l = const. 

Taking the k — > ±oo limit, we see that this constant has to be equal to 1, 
which finishes the proof. 

Proposition 4. The Ablowitz-Ladik scheme ( |1.17 ) with the parameters 



Oio = do = ct+ = S + = a_ = 0, 5- = 1 

is equivalent to the following map: 

(Qk - Qk)/h = q k+ i(l - q k r k ), 
T-\-h) : { (3.13) 
(r k - r k )/h = -r fc _i(l - q k r k ) 
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approximating the flow T+. This map has the commutation representation 



T-\-h) : L k wl +) = W { + + \L k 



with the matrix 

Wj k (+) =wi +) (q,r,/i) 




1 + hq k r k _i 

\ ""»' ft— 1 J- I '">iB' ft— 1 

'(3.14) 

Proof of this Proposition is omitted, because it is completely analogous 
to that of the previous one. 

Remark. It is important to notice the following relations: 



Wt\q,r,h) = (l-h\ 2 )(v k i+ \q,r,-h) 



WFfa&h) = (l + h\-'){Vr(cL,i,-h) / 

These relations are very difficult to guess and to prove, if one remains by the 
original formulation of the Proposition 0. 

We have demonstrated so far that the maps T + (h), T-(h) have com- 
mutation representations with the same matrix L k as their continuous time 
counterparts, and hence they share all the integrals of motion. To claim that 
the maps belong to the Ablowitz-Ladik hierarchy, we still need to show the 
Poisson property. 

Proposition 5. The both maps T + (h), T_(h) are Poisson with respect to 
the Poisson bracket ( |2.1|) . 

Proof. The Poisson property of a map (q, r) i— > (q, r) with respect to the 
bracket ( |2.1| ) is equivalent to the preservation of the corresponding 2-form: 

y, dq k A alr k _ " dq k A dr k 
^ 1 - q k r k ^ 1 - q k r k 

(remaining for simplicity by the finite dimensional case with the periodic 
boundary conditions). We shall prove this identity only for the map T + (h), 
since for T-(h) everything is completely analogous. Differentiating the equa- 
tions of motion (|3.1|), we obtain the following expressions: 



dq k = (1 - hq k+1 r k ) dq k + h(l - q k r k ) dq k+1 - hq k+1 q k dr k , (3.16) 
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(1 - hq h r k -i) dr k = dr k - h(l - q k r k ) dr k _ x + hr k _ x r k dq k . (3.17) 
Using succesively these two formulas, we obtain: 

dq k A dr k = (1 - hq k+1 r k ) dq k A dr k + h(l - q k r k ) dq k+ i A dr k 
1 - hq k+1 r k 



= ~ 7 — ~ — dq k A [dr k - h(l - q k r k ) dr k -i) + h(l - q k r k ) dq k+i A dr k . 

1 - hq k r k _i V ' 

(3.18) 

To make the last step, we observe that the equations of motion ( [3.1|) may be 
equivalently re-written as 

1 - qkh = (1 - qkh)(l - hqk+ih), (3.19) 

1 - qkTk = (1 - qkh){l - hq k h-i)- (3.20) 
Upon use of (|3~T9|) , (|33o| ) the identity ( gl| ) may be presented as 



dq k A dr fc _ A dr k h dq k+1 A dr k _ hdq k A dr fc _i 
1 - q k r k 1 - q k r k 1 - hq k+1 r k 1 - hq k r k -i ' 

Summation over k results in ( P-15| ). The proof is complete. 

It should be mentioned that each of the pairs (T + (h), Tf 1 (— h)) and 
(7^ 1 (— h), T-(h)) may be seen as generated by one of the two simplest par- 
titioned Runge-Kutta methods [[| when applied to the pair of differential 
systems (JF + ,JF_). Recall that the both of these partitioned Runge-Kutta 
methods are symplectic when applied to canonical Hamiltonian systems 



The Proposition 5 shows that this is still true for some systems which are 
Hamiltonian with respect to nonlinear Poisson brackets. However, unlike 
the canonical case, this statement cannot be extended to arbitrary systems 
Hamiltonian with respect to the bracket ( |2.1| ): the concrete form of the 
right-hand side is essential for the validity of the Proposition 5. 

We finish this Section by constructing a discretization for the linear flow 
JF . This is a much more simple task. Among many reasonable discretizations 
of the flow J-q we choose (for the reasons which will become clear in the next 
Section) 

%(h) : qk = ——rqk, r k = \^-rr k . (3.21) 
1 + h 1 — a 

(Note that T Q x {-h) = %(h)). 
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Proposition 6. The linear map % is Poisson with respect to the bracket 
J]) and has the commutation representation 



with the constant matrix 




(3.22) 



4 Local discretizations for T§ o T- o J r + 

Recall that, when considering the system fll.l| ) as a space discretization of 



the nonlinear Schrodinger equation ( |1.2|) , the following reduction is of the 
primary interest: 

r = ±q*, (4.1) 

(it is admissible in the case of pure imaginary values of time t, that is, after 
the change of the independent variable 1 1— > it, i — \J— 1). 

The flows JF + , JF_ alone do not allow this reduction any more, as well as 
their time discretizations T + (h), 71 (h). Nevertheless, we shall demonstrate 
now that the composition of these maps does again have this attractive prop- 
erty. 



Proposition 7. The Ablowitz-Ladik scheme ( |1.17|) with the parameters 
oco — So = ct- = 5- = 0, a + = 5 + = 1 
may be presented as the composition 

T-(h)oT + (h). 

Proposition 8. The Ablowitz-Ladik scheme ( |1.17|) with the parameters 
«o = o~o = ci + = 5 + = 0, a_ = 5_ = 1 
may be presented as the composition 

T^-h)- 1 oT + (-h)-\ 
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Proof of the Proposition 7. Let 

T+{h): (q,r)^(q,r), T.[h) : (q, r) h- (q, r), 
so that, according to fl3.1|), ( |3.6|) , 

(Qk - qk)/h = q k+1 (l - q k r k ), (r k - r k )/h = -f fc _i(l - g fe f fc ), (4.2) 

{qk-qk)/h = q k -i(l-q k r k ), (r k - r k )/h = -f k+1 (l - q k r k ). (4.3) 

It follows from the Propositions 1,2 that the composition 7L(h)oT + (h) allows 
the commutation representation 

L k V k = V k+l L k (4.4) 

with the matrix 

V k = Vt\^v)Vi + \^v). (4.5) 

We calculate now the entries of the matrix V k (denoting them according to 
( |1.12| )) in order to show that they have the form (|1.13|) — fll-16|) . From (|4.5|) , 



Q , (U we obtain: 

C k = -h\~ l r k {\ + h\ 2 - hq k r k -i) + h\f k (l - h\~ 2 + hq k -ir k ) 
= h\(r k -i - hr k (l - q k -ir k -i)) - h\~ l {r k + hr k -i(l ~ Qkh)) 
= hXrk-i - /iA _1 r fe 

(the last equality follows from the equations of motion (J4.2|), (J4.3|)), 

A k = 1 + h\ 2 - hA k 
V k = l-hX~ 2 + hD k , 

where 

A k = {q k + /ig fe _i)f fe _i, £> fc = - %)f fc . (4.6) 

So, the entries of the matrix V k has exactly the form ( |1.13| )-( "1.16|) with the 
parameters ao — do — a_ = 6- — 0, a + — S + — 1. We may conclude that 
the quantities A k , D k satisfy the difference relations ( |1.18|) , ( |1.19| ). (One 
could as well derive these difference relations (|1.18|) , ( |1.19|) directly from the 
definitions (|4.6|) and the equations of motion (|4.2|), (|4.3|).) In the case of 
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the rapidly decaying boundary conditions we have obviously A k} D k — > as 
k ±00, so that the quantities Ak, D k may be alternatively represented as 
in flOl]), ( [T22D . This finishes the proof. 

Proof of the Proposition 8. This time let 

TrH-h): fq,r)^(q,r), TT 1 ) 



-/i) : (q,r) i-> (q,r), 
so that, according to ( |3.13|) , (p^8|), 

{q k - q k )/h = q k+1 {l-q k r k ), {rk-r k )/h= -r k -i{l-q k r k ), (4.7) 

{qk-qk)/h = q k -i{l- qkr k ), {h - f k )/h = -r fc+ i(l - q k r k ). (4.8) 
It follows from the Propositions 3,4 that the composition Tz 1 {—h)oT^ l {—h) 
allows the commutation representation 

L k W k = W k+1 L k (4.9) 

with the matrix 



W k = Wt\q,r)W^'(q,r). (4.10) 

We calculate now the entries of the matrix W k (denoting them again accord- 
ing to ([L. 12 )) in order to show that they, as in the previous case, have the 
form ( [L.13 )-( 1.16 ). Denoting 

1 

A fc 



(1 - hq k -xr k )(l + hq k r k _ x 
we obtain from (WM , fl3jj), §] 
CfcAfc 1 = h\r k ^i - hX^r k 



(4.11) 



hXq k (l + hX~ - hq k -ir k ) - hX~ q k -i{l - hX + hq k r k _i) 
hX(q k + hq k -i{l - q k f k ] 
hXq k - hX~ l q k ^i 



hX - hq k (l - g fc _ir fc „i 



(the last equality being based on the equations of motion (|4.7|) , (|4.8|) ), 

1 - hq k _ x r k - h 2 q k _ x r k _ x + hX~ 2 



M 1 



Afc 1 - hr k ^i(q k + faf fc _i(l - g fc r fc )) + hX 



V h A 



Afe 1 - hq k r k _i + /iA 2 , 

1 + hq k r k -i - h 2 q k r k - hX 2 

Afc 1 + hr k (q k -i - hq k {\ - q k ^ x r k ^ 

A fc 1 + hq k _ x r k - hX 2 



-hX 2 
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(using again repeatedly the equations of motion ([4.7|), ( }4.8| )) and the defini- 
tion (PI|) ). 

We see that this time the entries of the matrix Wk have exactly the form 
( |1.13|) -( |I.16| ) with the parameters a = 5 = a + = S + = 0, a_ = 5_ = 1. 



It follows that the quantity A& satisfies the difference relation ( 1.20Q (which 



could be as well derived from the definition ( |4.11|) and the equations of motion 



(|4.7|) , (|4.8| ).) In the case of the rapidly decaying boundary conditions we have 
obviously — > 1 as k — > ±oo, so that may be alternatively represented 
as in (|1.23|) . The proof is finished. 

The maps constructed in these two Propositions have the desired property 
( |1.24|) , assuring that the reduction ( f4.1|) is admissible for them. However, 
they still do not approximate the Ablowitz-Ladik system ( |1 . 1| ) . In order to 
achieve this, we have to commute them with T (h). The following lemma, 
following directly from the formulas ( |1.17[ ), allows to control the parameters 
«o, Sq of the Ablowitz-Ladik discretizations. 

Lemma. Let T '(h) = 7al(^; 0, a+, a_; 0, 5 + , 5-), and let T"(h) be the 
linear map 

(q r ) i_> ( 1 _ q 1 + r \ 

\ 1 + h5o ' 1 — hao J 

T'(h) o T"(h) = T"{h) o T\h) = 

= T AL (h; a , a + (l - hao), - ha ); S , 8+(l + hS ), + hS )). 

According to this lemma, we derive from the Propositions 7,8 the follow- 
ing fundamental statement. 

Proposition 9. The Ablowitz-Ladik scheme ( |1.17 ) with the parameters 



Then 



cxq = Sq = 1, a + = 1 — h, 5 + = 1 + h, a_ = 5- = 
may be presented as the composition 

%{h)oT-(h)oT+(h). 
The Ablowitz-Ladik scheme ( |1.17| ) with the parameters 

oco = So = 1, a + = S + = 0, a- = 1 — h, 5_ = 1 + h 
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may be presented as the composition 



%(h) o Tz x {-h) o T+\-h). 



These two maps, approximating the system ( |1.17 ), have the property 



(|1.24|) under the condition h* = —h, i.e. h pure imaginary. Hence they 



both may serve as honest time discretizations of the reduced version (|1.5| ) 
of the Ablowitz-Ladik system ( |1.1| ). Note that all maps in each of these 
compositions commute. 



5 Connection with relativistic Toda 

As noted in [|ll|], the Ablowitz-Ladik hierarchy is in principle nothing but 
the relativistic Toda hierarchy and vice versa. We first give a Hamiltonian 
interpretation of this statement, and then use the results on the discrete 
time relativistic Toda lattice to clarify the place of the maps T + , 71 in 
the Ablowitz-Ladik hierarchy. 

Define new variables Ck, dk on the phase space of the Ablowitz-Ladik 
hierarchy: 

dk = , c fc = [q k r k - 1). (5.1) 

qk 

A direct computation shows that the only nonvanishing Poisson brackets 
between these functions are: 

{c k ,d k+ i} = Ckdk+i, {c k ,dk} = —Ckd k , {ck,Ck+i} = c k c k+1 . (5.2) 

One immediately recognizes in these relations the quadratic Poisson brackets 
underlying the relativistic Toda hierarchy. Moreover, one sees immediately 
that the simplest Hamiltonians of the Ablowitz-Ladik hierarchy H±, Hq may 
be expressed in the variables c k , d k as 

N N Ck _|_ d 

H + = J2 1k+ir k = a A * \ ( 5 ' 3 ) 

k=l k=l a fc«fc+l 

N N 

H- = qkTk+i = XX c <t + <^)> ( 5 - 4 ) 

k=l k=l 
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N 



N 



H = J2 log(l - q k r k ) = log(c fe /4), (5.5) 

k=l k=l 

and in H± we recognize the two basic Hamiltonians of the relativistic Toda 
hierarchy. We demonstrate now that our maps 7+, 71, being expressed in the 
variables Cfc, d k , also coincide with the discrete time flows of the relativistic 
Toda lattice introduced in 0. 

Proposition 10. Consider the map T + {h). Let the variables c k) d k be 
defined by (|5.1| ), and define the auxiliary function z>k by 

5fc = Qkh ~ 1- 
T/ien t/ie following relations hold: 

dk — h — fa>fc_i, 
4 — ^ffc-i ~ 



(5.6) 



(5.7) 



dk — dk+v 



Ck — Ck+V 



c k + hd h 



(5.8) 



dk+i — hdk Ck+i + hdk+i 

Proposition 11. Consider the map T_(h). Let the variables c k , dk be 



defined by (|5.l), and define the auxiliary function a k by 

Qk-i . hqk-i 



a-k 



q k -i 



(5.9) 



Then the following relations hold: 



1 + hd k + 



he, 



k-i 



dk = d, 



a-k+i 



- hd 



k+l 



Ck = c fc - 



Qfc+i + hck+i 



(5.10) 



(5.11) 



afc — hdk Ofe + hck 

Proof of the Proposition 10. We shall use in the proof the equations of 
motion in the form ( |3.19|) , ( |3.20| ) and also two additional auxiliary identities. 
The first equation of motion in ( |3 . 1| ) , re- written with the help of definitions 
for dk+i and dk, reads: 

— — = 4+i - hi k , (5.12) 

Qk+l 
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Further, from ( |3.19| ), ( |3.20| ) and the definition ( |5.6| ) it follows: 

1 - qkh _ <^k 

1 ~~ Qk+l r k+l f fc+l 



(5.13) 



Using the definitions ( |5.1| ) for Ck and ( p.6[ ) for then the identity ( |3.2(J| ), 
and again the definitions ( |5.1|) for dk and (|5.6| ) for Q k -i, we get: 



Cfc qk-i 1 — ?aA; <?fc-i / -i t — \ 7 j i 

— = — = (1 - hq k r k -i) =d k -h- hv k -i, 

*k qk 1 - gfc^fc gfc 

which is the recurrent relation ( |5.7| ). 

To prove the first equality in (|5.8|), we use the definition of d k and the 
identity (|5TT2|) : 

<4 _ qk-i qk+i _ d k — hdk-i 
dk+l qk qk dk+l — hD k 

Finally, to prove the second equality in (|5.8| ), we use the definition of c k , 
the identities ( |5.12| ), G5.13Q , and the recurrent relation (|5.7|) : 

c k qk-i q k +i 1 - Wk d k - /if fc _i 5 fc c k + hx> k 



Cfe+i qk qk 1 — qk+i^k+i dk+i — hd k dk+i Ck+i + hd 



k+l 



The Proposition 10 is proved. 

Proof of the Proposition 11. We start again with re- writing the 
equations of motion (|3.6| ) in the equivalent form: 



1 - q k r k = (1 - ?fc^)(l + hq k -ir k ), (5.14) 

1 - Mk = (1 - qkr k ){l + hq k r k +i). (5.15) 

Note that the first equation of motion in (|3.6| ) may be represented also in 
another equivalent form: 

q k (l + hq k -ir k ) = q k + hq k -i. (5.16) 



We shall need also two additional auxiliary identities. The definitions 
( |5.1| ) immediately imply: 

a k - hd k = < ^-, (5.17) 
qk-i 



19 



ak + hc k = -^-^{1 + hq k -ir k ) 
Qk-l 



(5.18) 



Now to prove the recurrent relation (|5.10|) we use (|5.18|) in conjunction 
with (pD, and then (|5l6|) and ( ^A7\) : 



1 + 



hc k q k (l + hq k _ir k ) 



a k 



qk + hq, 



k-l 



qk 

Qk 



hd 



fc+i- 



The first equation of motion in ( |5.11| ) follows from the definition of d k and 
the formula ( |5.17| ): 

dk _ Qk-i Qk _ gfc+i — hdk+i 
dk Qk-i q k a k — hd k 

Finally, to prove the second equation of motion in ( |5.11| ), we use the definition 
of c k , the formulas (|5.14|) , (|5.15| ), and then ( |5.16|) : 

Ck _ gfc-i Qk 1 - qkh _ gfc-i qk 1 + hq k r k+1 _ a k+1 + hc k+1 
Ck qk-i qk^~ qk^k ?fc-i <7fc 1 + hq k - X r k a k + hc k 

The Proposition 11 is proved. 

Now one immediately recognizes in ( |5.7| ), ( |5.8|) and in ( |5.10| ), ( |5.11|) two 
discrete time flows of the relativistic Toda hierarchy introduced and studied 
in 0. Applying the results in [Q], we get the following statement (formulated, 
as before, for the periodic case for the sake of notational simplicity). 

Proposition 12. The maps T± have Lax representations with either of 
the N x N Lax matrices 



T + (q,r,A) = L(q,r,A)t/- 1 (q,r,A) or T_(q, r, A) = U (q, r, A)L(q, r, A) 



where 



N _ N 

qk^i 



r ; A) = ^2 E kk + A ^ Ek+l,k, 

k=l Ik k=l 



A' 



N 



U(q,r, \) = E kk + \ 1 ^ -^(1 - q k r k )E k)k+ i. 

k=l k=l 1 k 

These maps are interpolated by the flows with the Hamiltonian functions 



-tr„($(-T± *)) =H+ + 0(h) and tr ($(T±)) = H_ + 0(h), 
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respectively, where 

$(0 = h- 1 / ? log(l + hrjfi- = £ + 0(h). 
Jo 77 

The initial value problems for the maps T± may be solved in terms of the 
matrix factorization problem for the matrices 

I - hT ± \t = 0)) and (l + hT±(t = 0)) , 

respectively. 

(In the formulation above tro stands for the free term of the Laurent 
expansion for the trace; the detailed definition of the matrix factorization 
problem in a loop group mentioned in this Proposition, may be found in @). 



6 Conclusion 

Re-considering the Ablowitz-Ladik discretizations from the modern point 
of view, undertaken in this paper, turned out to be rather fruitful. We 
factorized a highly non-local scheme into the product of very simple (local) 
ones, each of them approximating a more simple and fundamental flow of the 
Ablowitz-Ladik hierarchy. These local schemes may be stuided exhaustively. 
In particular, we found in this paper the interpolating Hamiltonian flows 
for them, as well as the solution in terms of factorization problem in a loop 
group. We guess that also in the practical computations our variant of the 
difference scheme will exceed considerably the old one. It would be interesting 
and important to carry out the corresponding numerical experiments. 

It seems also promising to re-consider from this point of view other non- 
local integrable discretizations derived and tested in 0. 

We note also that our maps are ideal building blocks for applying the 
Ruth-Yoshida-Suzuki techniques f[2" |, which will result in higher order in- 



tegrable discretizations for the Ablowitz-Ladik system. This point will be 
reported in detail elsewhere. 

The research of the author is financially supported by the DFG (Deutsche 
Forschungsgemeinshaft) . 
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